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Abstract 

A unified approach to (symmetric informationally complete) positive op- 
erator valued measures and mutually unbiased bases is developed in this arti- 
cle. The approach is based on the use of Racah unit tensors for the Wigner- 
Racah algebra of SU(2) D U(l). Emphasis is put on similarities and differ- 
ences between SIC-POVMs and MUBs. 
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1 INTRODUCTION 

The importance of finite-dimensional spaces for quantum mechanics is well recognized 
(see for instance |[Il-[l3]|). In particular, such spaces play a major role in quantum informa- 
tion theory, especially for quantum cryptography and quantum state tomography [|4ll- ll27l . 
Along this vein, a symmetric informationally complete (SIC) positive operator valued 
measure (POVM) is a set of operators acting on a finite Hilbert space [[4l|-[[T4| (see also 
|l3l for an infinite Hilbert space) and mutually unbiased bases (MUBs) are specific bases 
for such a space lfT5l-[|27ll. 

The introduction of POVMs goes back to the seventies ll4l|-|I7]|. The most general quan- 
tum measurement is represented by a POVM. In the present work, we will be interested 
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in SIC-POVMs, for which the statistics of the measurement allows the reconstruction of 
the quantum state. Moreover, those POVMs are endowed with an extra symmetry condi- 
tion (see definition in Sec. 2). The notion of MUBs (see definition in Sec. 3), implicit or 
explicit in the seminal works of [fT5l - [[T8l . has been the object of numerous mathematical 
and physical investigations during the last two decades in connection with the so-called 
complementary observables. Unfortunately, the question to know, for a given Hilbert 
space of finite dimension d, whether there exist SIC-POVMs and how many MUBs there 
exist has remained an open one. 

The aim of this note is to develop a unified approach to SIC-POVMs and MUBs based 
on a complex vector space of higher dimension, viz. (P instead of d. We then give a 
specific example of this approach grounded on the Wigner-Racah algebra of the chain 
SU(2) D U(l) recently used for a study of entanglement of rotationally invariant spin 
systems [|28l and for an angular momentum study of MUBs Il26ll27l . 

Most of the notations in this work are standard. Let us simply mention that I is the 
identity operator, the bar indicates complex conjugation, denotes the adjoint of the 
operator A, Sa^b stands for the Kronecker symbol for a and b, and A(a, 6, c) is 1 or 
according as a, b and c satisfy or not the triangular inequality. 

2 SIC-POVMs 

Let be the standard Hilbert space of dimension d endowed with its usual inner product 
denoted by ( | ). As is usual, we will identify a POVM with a nonorthogonal decompo- 
sition of the identity. Thus, a discrete SIC-POVM is a set {P^, : x = 1,2, ■ ■ ■ ,d'^} of d"^ 
nonnegative operators P^. acting on C^, such that: 

• they satisfy the trace or symmetry condition 

Ti {P,Py) = ^, x^y; (1) 

moreover, we will assume the operators Px are normalized, thus completing this 
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condition with 

Tr (P2) = 1; (2) 
they form a decomposition of the identity 

^E^^ = i; (3) 

x=l 

they satisfy a completeness condition: the knowledge of the probabilities px defined 
by Px = Tr(P^p) is sufficient to reconstruct the density matrix p. 



Now, let us develop each of the operators Px on an orthonormal (with respect to the 
Hilbert-Schmidt product) basis {ui : i = 1,2, ■ ■ ■ , cf} of the space of linear operators on 

Px = '^Vi{x)ui, (4) 

i=l 

where the operators Ui satisfy Tr (uluj) = 6ij. The operators Px are thus considered as 
vectors 

v{x) = {vi{x),V2{x), ■ ■ ■,Vd2{x)) (5) 

in the Hilbert space C^^ of dimension d"^ and the determination of the operators Px is 
equivalent to the determination of the components Vi{x) of v{x). In this language, the 
trace property ^ together with the normalization condition ^ give 

v{x)-v{y) = ^^idSx,y + l), (6) 

where v{x) ■ v{y) = Ylt=i '^i{x)vi{y) is the usual Hermitian product in C^^. 

In order to compare Eq. @ with what usually happens in the search for SIC-POVMs, 
we suppose from now on that the operators Px are rank-one operators. Therefore, by 
putting 

Px=\^x){^x\ (7) 
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with 10a;) e C^, the trace property ([Il[2l) reads 
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{dSx,y + 1) . 



d + 1 



(8) 



From this point of view, to find d'^ operators is equivalent to finding rf^ vectors 10^;) 
in satisfying Eq. ([8]). At the price of an increase in the number of components from 
d^ (for d"^ vectors in C'') to d"^ (for d"^ vectors in C^^), we have got rid of the square 
modulus to result in a single scalar product (compare Eqs. Q and ([8])), what may prove 
to be suitable for another way to search for SIC-POVMs. Moreover, our relation ^ is 
independent of any hypothesis on the rank of the operators P^. In fact, there exists a lot 
of relations among these d^ coefficients that decrease the effective number of coefficients 
to be found and give structural constraints on them. Those relations are highly sensitive 
to the choice of the basis {ui : i = 1,2, ■ ■ ■ , d"^} and we are going to exhibit an example 
of such a set of relations by choosing the basis to consist of Racah unit tensors. 

The cornerstone of this approach is to identify with a subspace e{j) of constant 
angular momentum j = {d — l)/2. Such a subspace is spanned by the set {|j, m) : m = 
—j, —j + 1, ■ ■ ■ , j}, where \j, m) is an eigenvector of the square and the 2;-component 
of a generalized angular momentum operator. Let u^'^) be the Racah unit tensor [|29ll 
of order k (with k = 0, 1, ■ ■ ■ , 2j) defined by its 2k + 1 components Ug'^^ (where q = 
—k, —k + l,---,k) through 



where (■ ■ ■) denotes a 3-jm Wigner symbol. For fixed j, the (2j + 1)^ operators Uq 
(with k = 0, 1, ■ ■ ■ , 2j and q = —k, —k + 1, ■ ■ ■ , A;) act on e{j) ~ and form a basis 
of the Hilbert space of dimension = (2j ' + 1)^, the inner product in being the 
Hilbert-Schmidt product. The formulas (involving unit tensors, 3-jm and 6-j symbols) 
relevant for this work are given in Appendix (see also [|29l to ffJll ). We must remember 
that those Racah operators are not normalized to unity (see relation (1461) ). So this will 
generate an extra factor when defining Vi{x). 




(9) 
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Each operator Pj. can be developed as a linear combination of the operators Ug''\ 
Hence, we have 

2j k 

p^ = J2Yl ' (10) 

k=0 q=-k 

where the unknown expansion coefficients Ckq{x) are a priori complex numbers. The 
determination of the operators is thus equivalent to the determination of the coefficients 
Ckq{x), which are formally given by 



Cfe,(x) = (2fc + l)($,|4')|$,), (11) 

as can be seen by multiplying each member of Eq. (flOl) by the adjoint of Up^ and then 
using Eq. (|46l) of Appendix. 
By defining the vector 

vix) = {vi{x),V2{x), ■ ■ ■ , vn{x)), N = (2j + 1)^ (12) 



via 
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Vi{x) = ^=^=Ckq{x), i = k'^ + k + q + l, (13) 
V 2A; + 1 

the following properties and relations are obtained. 



The first component f i(x) of v{x) does not depend on x since 

coo(x) = -i^=j (14) 
V2j + 1 

for all X G {1,2, ■■ -,(27 + 1)2}. 

Proof: Take the trace of Eq. (flOl) and use Eq. (l48l) of Appendix. 

The components Vi{x) of v{x) satisfy the complex conjugation property described 
by 



Ckq{x) = {-iyCk-q{x) (15) 

for all a; G {1, 2, ■ ■ ■ , (2j + 1)^}, A; g {0, 1, ■ ■ ■ , 2j} and q G {-k, -fc + 1, ■ ■ ■ , A;}. 
Proof: Use the Hermitian property of and Eq. (|43l) of Appendix. 



In terms of Ckq, Eq. © reads 

2j k 

J2 ^ ^^Ck,{y) = — ^ [(2j + 1)4,, + 1] (16) 

fe=0 q=-k ' 

for all x,y G {1, 2, ■ ■ ■ , (2j + 1)^}, where the sum over q is S0(3) rotationally 
invariant. 

Proof: The proof is trivial. 

The coefficients Ckq{x) are solutions of the nonlinear system given by 

fc=0 e=0 q=-kp=-i \ ^ ' ^ 

^ I ? ^ ^ i Cfc<j(x)Qp(a;) (17) 



_ J J J _ 

for all X e {1, 2, ■ ■ ■ , (2j+l)2}, fsT g {0, 1, ■ ■ ■ , 2j} andQ G {-iT, -/sT+l, ■■■,K}. 

Proof: Consider = and use the coupling relation (ISTI) of Appendix involving 
a 3-jm and a 6-j Wigner symbols. 

As a corollary of the latter property, by taking K = Q and using Eqs. (|47l) and (ISO)) 
of Appendix, we get again the normalization relation \\v{x) p = v{x) ■ v{x) = 1. 

• All coefficients Ckq{x) are connected through the sum rule 

(2i+l)2 2j k / ■ h ■ \ 

E E E ^^^(^) -m ! j^' = + 1) (18) 

a;=l fc=0 (}=-fc ^ ^ / 

which turns out to be useful for global checking purposes. 

Proof: Take the jm-jm' matrix element of the resolution of the identity in terms of 
the operators Px/C^j + 1). 

3 MUBs 

A complete set of MUBs in the Hilbert space is a set of d(d + 1) vectors \aa) G C"' 
such that 

\{aa\bf3)\^ = So.,pSa,b + ^(1 - 5a,b), (19) 



where a = 0,1, - ■ ■ ,d and a = 0, 1, ■ ■ ■ , d — 1. The indices of type a refer to the bases 
and, for fixed a, the index a refers to one of the d vectors of the basis corresponding to a. 
We know that such a complete set exists if is a prime or the power of a prime (e.g., see 
II61-II21). 

The approach developed in Sec. 2 for SIC-POVMs can be applied to MUBs too. Let 
us suppose that it is possible to find d + 1 sets Sa (with a = 0,1, - ■ ■ ,d) of vectors in C^, 
each set Sa = {\aa) : a = 0,1, - ■ ■ ,d — 1} containing d vectors \aa) such that Eq. (fT9l) 
be satisfied. This amounts to finding d{d + 1) projection operators 

liaa = \aa){aa\ (20) 

satisfying the trace condition 

Tr ijlajlbp) = Sa,f3Sa,b + ^(^ ~ ^'^^b), (21) 

where the trace is taken on C^. Therefore, they also form a nonorthogonal decomposition 
of the identity 

d d~l 
a=0 a=0 

As in Sec. 2, we develop each operator Yiaa on an orthonormal basis with expansion 
coefficients Wj(aa). Thus we get vectors w{aa) in C^^ 

w{aQ.) = {wi{aa) , W2{aa) , ■ ■ ■ , ^^2(00;)) (23) 

such that 

w{aa) ■ w{h(3) = 5a,p6a,b + ^(1 - ^a,b) (24) 
for all a, 6 G {0, 1, ■ ■ ■ , and a, /3 G {0, 1, ■ ■ ■ , d - 1}. 

Now we draw the same relations as for POVMs by choosing the Racah operators to be 
our basis in C^'. We assume once again that the Hilbert space is realized by e{j) with 
j = {d — l)/2. Then, each operator Haa can be developed on the basis of the (2j + 1)^ 
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operators m^^^ as 

2j k 

^-» = J2Y1 dkq{aa)uf\ (25) 

fe=0 q=-k 

to be compared with Eq. (flOl) . The expansion coefficients are 



dkq{aa) = {2k + 1) {aa\uq^'' \aa) (26) 

for all a G {0, 1, ■ ■ ■ , 2j + 1}, a G {0, 1, ■ ■ ■ , 2j}, A; G {0, 1, ■ ■ ■ , 2j} and g G {-k, -k + 
1, ■ ■ ■ , A;}. For a and a fixed, the complex coefficients dkq{aa) define a vector 

w{aa) = {wi{aa),W2{aa), ■ ■ ■ ,wn{0'O:)) , N = {2j + 1^ (27) 

in the Hilbert space C^, the components of which are given by 

Wi{aa) = dkg{aa), i = k'^ + k + q + l. (28) 

V2A; + 1 

We are thus led to the following properties and relations. The proofs are similar to those 
in Sec. 2. 



First component Wi^ao) of w{aa): 



doo{aa) = (29) 
V2j + 1 



for all a G {0, 1, ■ ■ ■ , 2j + 1} and a G {0, 1, ■ ■ ■ , 2j}. 
Complex conjugation property: 



dkq{aa) = (-l)''4_g(aa) (30) 

for all a G {0, 1, ■ • ■ , 2j + 1}, a G {0, 1, ■ ■ ■ , 2j}, k G {0, 1, ■ ■ ■ , 2j} and q G 
{-k,-k + l,---,k}. 

Rotational invariance: 

5^ 2k + 1 ^ dkq{aa)dkq{h(3) = 5a,i35a,b + ^(1 " ^^.f*) (31) 

k=0 q=-k 

for all a, 6 G {0, 1, ■ ■ ■ , 2j + 1} and a, /3 G {0, 1, ■ ■ ■ , 2j}. 
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Tensor product formula: 



fc=0 ^=0 q=-kp= 

( k i K \ 

X I j j j J dkq{aa)dfp{aa) (32) 
for all a G {0, 1, ■ ■ ■ , 2j + 1}, a G {0, 1, ■ ■ ■ , 2j}, K G {0, 1, ■ ■ ■ , 2j} and Q G 

• 5Mm rule: 

2j+l 2j 2j k / ■ h ■ \ 

E E E E ^^^(««) Urn a m' = (-l)'""^2(2j + (33) 

a=0 a=0 fc=0 q=-k ^ ^ / 

which involves all coefficients dkq{aa). 

4 CONCLUSIONS 

Although the structure of the relations in Sec. 1 on the one hand and Sec. 2 on the other 
hand is very similar, there are deep differences between the two sets of results. The 
similarities are reminiscent of the fact that both MUBs and SIC-POVMs can be linked to 
finite affine planes [[l2l[Il|22l|2l|25l and to complex projective 2-designs lf8l[T0l[T9ll24l. 
On the other side, there are two arguments in favor of the differences between relations © 
and (l24l) . First, the problem of constructing SIC-POVMs in dimension d is not equivalent 
to the existence of an affine plane of order d [fT2l[T3l . Second, there is a consensus around 
the conjecture according to which there exists a complete set of MUBs in dimension d if 
and only if there exists an affine plane of order d [|22l . 

In dimension d, to find d'^ operators of a SIC-POVM acting on the Hilbert space 
amounts to find d'^ vectors v(x) in the Hilbert space with N = d'^ satisfying 

ll^^ll = 1, v{x) ■ v{y) = -j^-^ioT X ^ y (34) 

(the norm ||f of each vector v{x) is 1 and the angle tOxy of any pair of vectors v{x) 
and v{y) is tUxy = cos~^[l/{d + 1)] for x ^ y). 



In a similar way, to find d + I MUBs of C"' is equivalent to find d + 1 sets Sa (with 
a = 0,1, - ■ ■ ,d) of d vectors, i.e., d{d + 1) vectors in all, w(aa) in with N = d^ 
satisfying 

w^aa) ■ w{a(]) = Sa,i3, w{aa) ■ w{bf3) = - for a 7^ 6 (35) 

(each set Sa consists of d orthonormalized vectors and the angle tOaabp of any vector 
w{aa) of a set Sa with any vector w(6/5) of a set Sb is ujaabp = cos''^{l/d) for a 7^ b). 

According to a well accepted conjecture (H [TOl, SIC-POVMs should exist in any 
dimension. The present study shows that in order to prove this conjecture it is sufficient 
to prove that Eq. (|34|) admits solutions for any value of d. 

The situation is different for MUBs. In dimension d, it is known that there exist d + 1 
sets of d vectors of type \aa) in satisfying Eq. (fT9l ) when is a prime or the power of 
a prime. This shows that Eq. (|35l) can be solved for d prime or power of a prime. For d 
prime, it is possible to find an explicit solution of Eq. (fT9l) . In fact, we have [|26ll27l 

1 ^ 

\aa) = ^ J2 m), (36) 

* m=—j 

for a, a G {0, 1, ■ ■ ■ , 2j} while 

\aa) = \j, m) (38) 



for a = 2j + 1 and a = j + m = 0,1, ■ ■ ■ ,2j . Then, Eq. (1261) yields 

j j 



- Iri E -"'-''(-ir- ( 4„ J ^, ) . (39) 



m=—j m'=—j 

6{m,m') = {m — m') 



-(1 — m — m')a + a 



(40) 



for a, a G {0, 1, ■ ■ ■ , 2j} while 



d,,{aa) = 6,,o{2k + l){-iy-"^[ 4„ q ^ ) (41) 



for a = 2j + 1 and a = j + m = 0, 1, ■ ■ ■ , 2j. It can be shown that Eqs. (1401) and (1411) 
are in agreement with the results of Sec. 3. We thus have a solution of the equations for 
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the results of Sec. 3 when d is prime. As an open problem, it would be worthwhile to find 
an explicit solution for the coefficients dkq{aa) when = 2j + 1 is any positive power 
of a prime. Finally, note that to prove (or disprove) the conjecture according to which a 
complete set of MUBs in dimension d exists only if is a prime or the power of a prime 
is equivalent to prove (or disprove) that Eq. (|35l) has a solution only if ci is a prime or the 
power of a prime. 

APPENDIX: WIGNER-RACAH ALGEBRA OF SU(2) D 
U(l) 

We limit ourselves to those basic formulas for the Wigner-Racah algebra of the chain 
SU(2) D U(l) which are necessary to derive the results of this paper. The summations 
in this appendix have to be extended to the allowed values for the involved magnetic and 
angular momentum quantum numbers. 

The definition dH) of the components uf^^ of the Racah unit tensor u^^^ yields 

(j,m|«f |j,m') = (-l)^--™( 4^ g ^')' (42) 

from which we easily obtain the Hermitian conjugation property 

uf = i-iyu^^l (43) 

The 3-jm Wigner symbol in Eq. (I42l) satisfies the orthogonality relations 

Y ( ^ ^ V -^'z M = ^^4ApA(j,/,A;) (44) 

mm' ^ ^ ^ ^ 

and 

*)(m M' *)=W„..,m.. (45) 

The trace relation on the space e{j) 

Tr (y^^vf) = ^^4A,pA(j, J, k) (46) 
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easily follows by combining Eqs. (l42l) and (|44l) . Furthermore, by introducing 

in Eq. (I44l) . we obtain the sum rule 

E(-l)'"'" ( -"L g ^) =v/2jTT4AoA(j,A:,j), (48) 

known in spectroscopy as the barycenter theorem. 

There are several relations involving 3-jm and 6-j symbols. In particular, we have 



mm' M 



3 k j \ / j i j \ f j K j 
—m q M I \ —M p m' J \ — m Q m' 



where {• ■ ■} denotes a 6-j Wigner symbol (or W Racah coefficient). Note that the intro- 
duction of 

^ ^ ° \ 1 (50) 

i W J ' V(2A; + l)(2j + l) 

in Eq. (|49l ) gives back Eq. (|44l) . Equation (|49l ) is central in the derivation of the coupling 
relation 



ufu. 



Equation (|5TI) makes it possible to calculate the commutator [u'q\up''] which shows that 
the set {u"q^ : A; = 0, 1, ■ ■ ■ , 2j; q = —k, —k + 1, ■ ■ ■ , A;} can be used to span the Lie 
algebra of the unitary group U(2j + 1). The latter result is at the root of the expansions 
dn]) and (1321). 

Note added in version 3 

After the submission of the present paper for publication in Journal of Russian Laser 
Research, a pre-print dealing with the existence of SIC-POVMs was posted on arXiv [132| . 
The main result in [l32ll is that SIC-POVMs exist in all dimensions. As a corollary of this 
result, Eq. (34) admits solutions in any dimension. 
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